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In this paper we extend the work synthetically presented in Ref. [1] and give theoretical details 
I and complete tables of numerical results. We exploit calculations within a Bethe-Salpeter (BS) 

■ formalism adjusted for QCD, in order to extract an "experimental" strong coupling (^^^(Q^) below 
' 1 GeV by comparison with the meson spectrum. The BS potential follows from a proper ansatz on 

OA , the Wilson loop to encode confinement and is the sum of a one-gluon-exchange and a confinement 

^jQ' terms. Besides, the common perturbative strong coupling is replaced by the ghost-free expression 

r—{ ' aE(Q^) according to the prescription of Analytic Perturbation Theory (APT). 

^^rfH , The agreement of af'^{Q^) with the APT coupling a^iQ^) turns out to be reasonable from 1 GeV 

^ ' down to the 200 MeV scale, thus confirming quantitatively the validity of the APT prescription. Be- 

' low this scale, the experimental points could give a hint on the vanishing of aa((3^) as Q approaches 

zero. This infrared behaviour would be consistent with some lattice results and a "massive" gener- 
' alization of the APT approach. 

(— I , As a main result, we claim that the combined BS-APT theoretical scheme provides us with a rather 

satisfactory correlated understanding of very high and rather low energy phenomena from few hun- 
I ' dreds MeV to few hundreds GeV. 

^: 
I) 

^ . 

' ^' I. INTRODUCTION 

(N ■ 

This paper is an extended and more thorough analysis of the work already presented concisely in Rcf. [1]. Here 
^\ we try a self-contained theoretical discussion of the method, with all the premises and the necessary details. We 

■ report the complete numerical tables of our results, that were originally displayed only graphically, and we discuss 
'""^ I them extensively. As it is well known, a rather consistent picture of many strong high energy processes was obtained 

■ by perturbative QCD, if the running coupling a^iQ^) , as derived from the renormalization group equation, is used. 
_ Inversely, if the QCD scale value A„^=5 ~ 200 MeV is taken, the values of extracted from the data for the 
I appropriate Q or y/s fit rather well the theoretical as{Q^) curve, with few exceptions. Quite important, the 2- loop 

level of approximation for as seems to be sufficient [2] for practical size of data errors. 

Unhappily, as{Q^) develops (at any loop level) unphysical singularities for Q ~ ^nf=3 ~ 400 MeV, that makes the 
expression useless in the infrared region. This is a particularly serious difficulty in any type of potential model in 
^ , which Q should be identified with the momentum transfer, that takes typically values below 1 GeV, according to the 

■ state and to the mass of the quarks implied. 

Among various proposals to eliminate these singularities (see, e.g., Section 3 in Ref. [3] ), we mention here two particular 
ones, i.e., the freezing hypothesis, that simply consists in assuming that cXsiQ^) freezes to a certain maximal value H 
in the infrared (IR) region as a consequence of non-perturbative effects, and the analytization prescription of the APT 
approach [4]. The latter requires as{Q^) to satisfy a dispersion relation with the only unitary cut for — cx3 < < 
(throughout this paper we assume the momentum scale to be spacelike = —q^ > 0) and uses perturbation theory 
to evaluate the spectral function. 

On the other hand, in the last years various relativistic formalisms have been proposed in the context of QCD, or 
QCD motivated, that take confinement into account and evaluate the meson (and baryon) spectrum in the light and 
in the heavy quark sectors. Among the most recent works we remind, e.g., [5-10] and references therein. In this 
paper we reverse somewhat the point of view, that is, we take advantage of the comparison between the calculated 
meson spectrum and the data in order to gain information on the infrared behaviour of the QCD coupling that we 
shall compare with APT. To this end, we exploit a Bethe-Salpeter (BS) like formalism (second order BS formalism) 
developed in [11] and applied with a certain success to the calculation of a rather complete quarkonium (meson) 
spectrum in Refs. [12,13]. The formalism is essentially derived from the QCD Lagrangian taking advantage of a 
Feynman-Schwinger representation for the solution of the iterated Dirac equation in an external field. Confinement 
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is encoded through an ansatz on the Wilson loop correlator; indeed the expression i In is written as the sum of a 
one-gluon exchange (OGE) and an area term 

ilnW = {ilnW)oGE + (tS . (1) 

By means of a three dimensional reduction, the original BS equation takes the form of the eigenvalue equation for a 
squared bound state mass 

= Ml + [/oGE + ?7conf , (2) 

where Mq is the kinematic term Mq = w\+W2 = y'mf + + y^m| + k^, k being the cm. momentum of the quark, 
mi and the quark and the antiquark constituent masses, and U = Uqge + f^Conf the resulting potential. 

As a consequence of ansatz (1), the perturbative part of the potential C/qge turns out to be proportional to as{Q'^) , 
where in a sense as{Q^) should be identified as an effective charge of the type proposed in [14] and denoted in [3] 
with asGD(Q^)- 

Calculations have been performed in Refs. [12,13] by using both a frozen and the 1-loop analytic coupling (^^''(Q^) 
with an effective scaling constant A^^^^'^' ~ 200 MeV (see Eq. (9) below), which is equivalent at the 3-loop level to 
Ai%3 - 400 MeV or to the world average A^f^^^g ~ 200 MeV. 

The results of the two sets of calculations are relatively similar for the heavy-heavy quark states. However, for the 
IS states involving light and strange quarks, quite different results have been obtained in the two cases. In the 
case of a frozen coupling the tt and K masses turn out to be too high, independently of how small the light quark 
mass is taken (see Fig. 1); e.g., if we fit the light and the strange quark masses to the p and the (p masses, we find 
TOtt ^ 500 MeV and niK ^ 700 MeV, respectively. On the contrary, if appropriate values for the quark masses are 
chosen, the tt, p, K, K*, (j) masses can be rather well reproduced when the analytic coupling a^\Q'^) is used. This 
occurrence strongly supports the use of the analytic coupling in the BS framework. 

In this work we undertake a thorough analysis of this issue, by comparing our theoretical results, obtained for a certain 
choice of the parameters and the analytic coupling, with the results of a similar calculation performed by means of a 
fixed value of as, for every quark-antiquark state. We denote by af^ the value that reproduces the same theoretical 
result as obtained with oi"^\q'^) and by the corresponding value that reproduces the experimental mass. The 
value a^^ is then used to identify an effective Q pertaining to that particular state, which is to be understood as the 
argument of the related "experimental" coupling al^'^{Q'^). 

Since only the leading perturbative contribution in the BS kernel has been included, a rough estimate of NLO effects 
on the ctg^P value leads to a relative theoretical error which spans from 20% to much less than 1% throughout the 
spectrum according to the quark masses involved. Furthermore, since coupling among different quark-antiquark 
channels has not been taken into account, the theoretical masses are expected to reproduce the experimental ones 
within the half width r/2 of the state. In the framework of the BS formalism these are the most relevant sources of 
theoretical error and overwhelm all other errors, like those related to the three dimensional reduction. When relevant, 
the experimental error, related to the uncertainty of the experimental mass is added to the theoretical one. 
It should be noted that our results are model dependent in the sense that they do depend on the ansatz (1). More 
sophisticated models also exist, like the Stochastic vacuum model [5] and dual QCD [6], but they turn out to be 
considerably difficult to implement in the BS formalism. 

If compared with the 3-loop analytic curve a^\Q'^) normalized at the world average ag''(M|) = 0.1176 (A^^^''^^ = 
236 MeV), the values of af^^iQ"^) fit it rather well within error bars in all the region Anj^3/2 <Q < SA^^j^^ (being 
A^^^''^3 = 417 MeV according to threshold matching). In the region Q < our data fall below a^\Q'^) , and, 

with the limitation due to the large errors, this could be interpreted as a hint on the vanishing of q;s(Q^) as Q ^ , 
or on the existence of a finite limit lower than the universal Shirkov-Solovtsov freezing value (see e.g. [9]). Note, 
however, that, as it will be discussed extensively in Sec. 5, the experimental situation is particularly uncertain in this 
region and the theoretical treatment problematic. Nonetheless, it is worthwhile to mention that the former behaviour 
is consistent with a recently developed "massive" version of the analytic approach for the QCD coupling [15] (see Sec. 
2), and some results from lattice simulations. 

Finally let us stress that the choice to compare af^'P{Q'^) with the 3-loop expression a^\Q^) was to stay as close as 
possible to the usual practice in perturbation theory. In APT, however, when the appropriate small change of scale 
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(2) (3) 

is made (A^^'^g ~ 258 MeV rather than ^nf=5 — ^36 MeV if both normalized at the Z mass*) the 2-loop couphng 

ag^((5^) is practically indistinguishable from a^\Q^) in the entire Euclidean range and can be used instead of the 
latter for all practical purposes. 

The layout of the paper is as follows. In Sec. 2 the ghost-pole problem is discussed and an overview of the key 
points of analytic approach to QCD is given. In Sec. 3 an explicit expression for the BS potential U is given and 
the mathematical method to treat the eigenvalue problem for the squared mass operator is described. Sec. 4 is 
devoted to the strategy for extracting af^^iQ"^) from the data and to errors estimate. Finally, in Sec. 5 our results 
are discussed and a match of the extracted QCD coupling with relevant high energy experimental data is attempted 
via Analytic Perturbation Theory. 

Some technical material is exposed in Appendices. A brief review of the derivation of the second order BS formalism 

and of the expression of from the ansatz (1) is given in App. A. Numerical tables in App. B display all results 
in details. In App. C an useful formula for 3-loop analytic coupling (the spectral density) is explicitly given and 
compared with the usual 3-loop perturbative coupling. 



II. ANALYTIC APPROACH TO QCD 

The renormalization group (RG) method is an inherent part of theoretical description of strong interaction processes. It 
is usually employed to improve the results of perturbation theory in the high energy region. However, a straightforward 
application of the RG method to perturbative expansion eventually gives rise to unphysical singularities of both the 

RG-invariant coupling function^ and physical observables. The presence of these singularities contradicts the general 
principles of the local QFT^ and severely complicates theoretical analysis of hadron dynamics in the IR domain. 

There is a variety of the nonperturbative approaches to handle the strong interaction processes at low energies (for a 
recent review of this issue see Sec. 3 in paper [3] and references therein). Some of such methods originate in the general 
properties of perturbative power series for the QCD observables in the framework of the RG formalism. For instance, 
these are the method of effective charges [18], the "optimal conformal mapping" method [19] (see also Rcf. [20]), 
the "optimized perturbation theory" [21], and some others. There is also a number of approaches which impose 
nonperturbative constraints either on the strong running coupling (see, e.g., Refs. [22,23]) or on the RG (3 function 
(see, e.g., Refs. [24-26]). In this paper we will exploit the so-called Analytic Perturbation Theory (APT) [4] and its 
recent "massive" modification [15], which are briefly overviewed in what follows. 

The analytic approach to QFT constitutes the next step (after the RG-summation) in improving the perturbative 
results. Specifically, in addition to the property of renormalizability this method retains a general feature of local QFT, 
the property of causality. The analytic approach has first been devised in the context of Quantum Electrodynamics [27] , 
and then extended to the QCD case about ten years ago [4]. The basic merits of the analytic approach to QCD are 
the absence of unphysical singularities of the invariant charge and the enhanced stability of outcoming results with 
respect to both higher loop corrections and choice of renormalization scheme [28]. Besides, this method enables one 
to process the spacelike and timelike data in a congruent way [29]. A fresh review of the analytic approach to QCD 
and its applications can be found in paper [30] (a generalization of APT for fractional powers of ag was implemented 
in Ref. [31]). 

Usually, in the framework of RG-improved perturbation theory a QCD observable D{Q'^) of a single argument 
= —q^ >Q (spacelike momentum transfer squared) can be represented as power series in the strong coupling as{Q^): 

i?p,(Q2) = i + ^d„[a,(Q2)]", (3) 

ra>l 



*Note that these values of the scale constant turn out to be somewhat larger than the perturbative values, as given e.g., by 
Bethke [2] , with the same normalization for Qs . 

^For example, the one-loop strong running coupling (4) possesses the so-called Landau (or spurious) pole in the low energy 
region. This problem can not be solved by the inclusion of higher loop corrections since the latter just give rise to additional 

singularities of the cut type. 

*It is worth noting also that the results of lattice simulation testify to the absence of spurious singularities of the QCD couphng 
at low energies, see, e.g., a recent overview in Sec. 2 of paper [16] as well as original papers [17]. 
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where dn arc the relevant perturbative coefficients. However, in the IR domain this expansion becomes inapplicable 
due to spurious singularities of the running coupling as{Q^). For example, the one-loop expression 

possesses both the physical cut along the negative real semiaxis <Q and unphysical pole at = A^. 

In the framework of the APT, the power series (3) for an "Euclidean" observable is replaced [32] by the nonpower 
expansion 

D^,,{Q^) = l + J2dnAn{Q^) (5) 

n>l 

over the set of functions 

CXD 





Here, the spectral function pia) is defined as the discontinuity of the relevant power of the perturbative running 
coupling ctsiQ^) across the physical cut, namely 

Pn(o-) = - Im[as{-a - is)]"' . (7) 

TT 



The APT representation for a QCD observable D{Q'^) (5) is free of spurious singularities. Besides, it displays a better 
stability (in comparison with the perturbative parameterization (3)) with respect to both, higher loop corrections and 
choice of the renormalization scheme, see Ref. [30] for details. The first-order function Ai{Q^) (6) plays the role of 
the effective Euclidean QCD coupling at respective loop level: 

oo 

aE(O') ^ = J da. (8) 





In the one-loop case this equation can be integrated explicitly [4] 

1 



Po 



A2 



ln(QVA2) A2-g2 



(9) 



At the higher loop levels the spectral functions (7) become rather involved (see App. C). An extensive numerical 

study of the analytic running coupling (8) and its "effective powers" (6) at various loop levels can be found in 
Ref. [33]. Besides, for practical applications one can also use simple explicit expressions [34] which approximate the 
APT functions (6) within reasonable accuracy. 

In order to handle the QCD observables which do not satisfy the integral representation of the form of Eq. (6), the 
APT has to be modified appropriately. For example, the Adler function, being defined as the logarithmic derivative 
of the hadronic vacuum polarization function, satisfies the dispersion relation [35] 

oo 
4ml 



where R{s) denotes the Drell ratio of the electron-positron annihilation into hadrons. Thus, the Adler function (10) 
can be expanded over the set of functions AniQ"^) (6) only in the limit of the massless pion m„ = 0, since otherwise 
the analytic properties of D{Q'^) (10) in variable differ from those of An{Q'^) (6). 
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The effects due to the nonvanishing pion mass have been incorporated into the analytic approach to QCD in Ref. [15]. 
In the framework of the latter formalism the Adler function (10) can be expanded 

D^,,^{Q\ml) = -^^^^ + J2dnK{Q\ml) (11) 



over the set of the "massive" functions 

a - Am? da 

4m^ 



Obviously, in the massless limit Eqs. (11) and (12) coincide with the expressions (5) and (6), respectively. Similarly 
to the case of the APT (8), the first-order function Ai(Q^,m^) (12) plays the role of an effective "massive" running 
coupling at the relevant loop level, namely 

oo 

"W^-') -e^^;? / (13) 

imp- 

It is worthwhile to note that irrespective of the loop level this coupling possesses the universal IR limiting value 
a.{Q'^,-rr?) ^ at = 0, see Ref. [15] for details. 



III. BS-MODEL FOR QQ SPECTRUM 



As mentioned, in [12,13] the meson spectrum is obtained by solving the eigenvalue equation for the squared mass 
operator (2), where the perturbative and confinement part of the potential are respectively 

(k|C/oGE|k') = 



4as(Q^) /(wi +W2)(wi +w^) 



3 7r2 

i 



^ 2 (Q ■ q)' 



2 k X k' • (ai + 0-2) 



2Q 

1 1 (\ 

+ g^i-a2 + -(-ar-a2 



920 ("1 • Q) - 910 ("2 • Q)] + 

2(ai-Q)(a2-Q) 



2Q2 

(Q-^Tl)(Q-fT2)\ , 1 



Q2 



4Q2 



(14) 



and 



(27r)3 



j d3re''i '-J"-Hi-,q,gio,«2o) 



(15) 



with 



J'""'(r,q, 910,920) 



910 + 920 



920 V 9?o - qi + 9?o V «io - qi 



^2 „2 



910920 / . q± , . q± 

- arcsm h arcsm 

|q±l V 910 920 



920 



+ 



V^ito - q± 

910 



(r X q - fji +i9io(r • ai)) 



- qi 



(r X q - f72 - i92o(r • 02)) 



(16) 
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Here denote the usual Dirac matrices 7j7j, cr^ the 4x4 Pauli matrices and q = , Q = 

k — k' , Qjo = — '^1 "^2 are constituent masses. 

Eqs. (14-16) follow from the ansatz (1) and a 3-dimcnsional reduction of our Bctho-Salpotcr equation (see Refs. [11-13] 
and App. A for the details). Actually in the calculation of [12,13] only the center of gravity (c.o.g.) masses of the fine 
multiplets were considered as a rule, and the spin dependent terms in (14-16) (spin-orbit and tensorial terms) were 
neglected with the exception of the hyperfine separation term in (14) proportional to |(Ti • (72 . Within this limitation 
a generally good reproduction of the spectrum was obtained for appropriate values of the parameters, as apparent 
from Fig. 1 re-elaborated from [13]. Here the results of three sets of calculations are displayed. Diamonds refer to 
the usual perturbative 1-loop coupling (to be replaced in Eq. (14)), frozen at a maximum value H , which has been 
taken as an additional adjustable parameter. Squares and circles both refer to the 1-loop APT coupling (9) with 
A ~ 200 MeV. For light quarks a running constituent mass was used too.^ 

We stress that only with the choice (9) the l^^o state has been correctly reproduced when light and strange quarks 
were involved, as in the case of tt and the K mesons. 

In the present work a similar calculation with the input (9) and a slight different choice of the parameters is made 
(preliminary results were given in [36]). First, we have fixed the string tension to the value cr = 0.18 GeV^ (consistent 
with other phenomenology and lattice simulations) and the scale constant to a!^''^J^] = 193 MeV. The whole set of 
remaining parameters, all the quark masses, are then determined by fitting the tt , (p , J/ip and T masses. It turns 
out m„ = nid = 196 MeV, = 352 MeV, = 1.516 GeV and m;, = 4.854 GeV. The results for the meson spectrum 
are given in the fourth column of tables in App. B. 

The chosen value for A^^^'^^ has been dictated by the comparison with the 3-loop analytic coupling normalized at 
the Z boson mass (see App. C) according to the world average. As displayed in Fig. 2 the relative difference between 
the two curves is no more than 1% in the region 0.5 < Q < 1.2 GeV, to which the states used as an input in the 

calculation belong. 

Furthermore, as already noted, our equations refer to a single definite quark-antiquark channels. So, having correct 
relativistic kinematics, they do not include coupling with other channels like any potential model (see App. A). Then 
we can not expect to have any insight into the splitting of over-threshold complicated multiplets which involve mixture 
of different states. Even the position of the c.o.g. mass is expected to be reproduced only within one-half of the width 
of the state. This has been taken into account in the estimate of the theoretical error (see Sec. 4). 
The resolution method of the eigenvalue equation for the operator (2, 14-16) we have used in [12,13] and in the present 
work can be summarized in the following way. 

a) In the static limit the problem can be reduced to the corresponding one for the center of mass Hamiltonian (see 
App. A) 

4 (X 

HcM = wi +W2 - - — + ar . (17) 



b) The eigenvalue equation for (17) is solved for a convenient fixed as by the Rayleigh-Ritz method, using the three 
dimensional harmonic oscillator basis and diagonalizing a 30 x 30 matrix. 

c) The square of the meson mass is evaluated as ((/)o]M^|(/)a), 4'a being the eigenfunction obtained in step b) (with a 
the whole set of quantum numbers) and the operator given by Eq. (2). 

d) Prescription c) is equivalent to treat — as a first order perturbation. Consistently the hyperfine separation 
should be given by 

Crrinif - Crrinif = — dkk^ dk' k'\li{k)^^i{k') x 



^Circles refer to a phenomenological running mass as function of the cm. quark momentum = = 0.17|k| — 0.025|k|^ -|- 
O.lSjkj'*. Squares refer to a running constituent mass resulting from a solution of the DS equation with an analytic RG running 
current mass (sec App. A for the details) which was, however, more in line with an attempt to define an analytic a^\Q^) 
singular at — > and including confinement [24,25] than to (9). 
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IWi + W2 W[+W2 

W1W2 Y W[w'2 J_i 



d^a^mPiiO, (18) 



where (pni is the radial part of the complete eigenfunction (j)a ■ However, in the case of the states involving light and 

strange quarks the quantity is further corrected to the second order of perturbation theory. 

For the quark masses and string tension a in (17) we have used the same values listed above and for what regards ag , 
that is supposed to be a constant in (17), we have taken as = 0.35 , which is the typical value used in non-relativistic 
calculations and also the freezing value adopted in [12]. 



IV. EXTRACTING ai^^{Q'^) FROM THE DATA 



One focus now on the reversed problem, i.e., the determination of the af^^{Q'^) values at the characteristic scales of 
a selected number of ground and excited states. 

In order to estimate af^^{Q'^) at low scales one needs first to assign an cfFcctivc Q-value to each state. To this end 
one first rewrites the squared mass, as given by point c) in Sec. 3, more explicitly as the sum of the unperturbed part, 
the perturbative and the confinement one respectively 

ml = {<Pa\M^\^a) + ('/'a|C/oGE|0a) + (0a |C/Conf ■ (19) 



Here J/qge is given by the second line of (14) and J/conf by Eq. (15) and first two lines of (16). Prom the OGE 
contribution we then extract for each state the fixed coupling value a^^ which leads to the same theoretical mass as 
by using Qg^(Q^) given by Eq. (9). This can be done by means of the relation 

(<^a|t^OGE|.^a) = (</.„|4'^ (Q2)0(q; Q)|<^„) = af (</.„|0(q;Q)|.^„), (20) 



where ^^(q; Q) can be drawn again by the second line of Eq. (14). The effective momentum transfer Qa associated to 
each bound state is then identified by equating 

4'^(qD = "1^ (21) 

The next step is to search for the correct (fixed) value of the coupling that exactly reproduces the experimental mass 
of each state. This is defined by the relation 

{cl>a\M^\cl>a) + at-PiQl){cl>M<l^ QMa) + {ct>a\Ucon{\<l>a) = ml^^ , (22) 



so that, by combining Eqs. (19), (20) and (22) we finally obtain 

1 0(q;Q) 



^exp ) ^ ml^^-ml + af{^a\OM)\cl>a) _ ^^^^ 



This procedure has been applied to a number of light-light, light-heavy and heavy-heavy ground as well as excited 
states. 

Note that, apart from the particular ansatz made in (1) to take into account confinement, the other relevant approx- 
imations are: 



i) only the leading perturbative contribution is included in (1) and so in the potential; 

ii) quark antiquark annihilations and couplings with other channels have been ignored; 

iii) an instantaneous approximation is involved in deriving the eigenvalue equation for (2) from the original BS equa- 
tion. 



As even the experience with QED suggests, retardation corrections are expected to be relevant for the hyperfine 
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and possibly fine splitting, but of minor importance for the positions of the c.o.g. of the multiplets that we essentially 

use to evaluate q!™p((5^). 

Thus, as we told, the main sources of theoretical error in the whole procedure arc expected to arise from neglecting the 
NLO contribution to the BS kernel as well as the coupling with other channels. For what concerns the former (point 
i)), it is worth noting that the next to leading contribution to the perturbative part of the BS-kernel comes from four 
diagrams with two-gluon exchange; two triangular graphs containing a four- line vertex of the type g'^cj)* cfiA^A'^ and 
two three- line vertices g<l)*d^<pA'^ (the spin independent part of our second order formalism is quite similar to scalar 
QED), one fish diagram with two four-line vertices, and a crossing box with four three line vertices. If the renormal- 
ization scale is identified with the momentum transfer Q the fish graphs contribution is completely reabsorbed in the 
renormalization. On the other hand, a somewhat crude estimate of the contribution of each of the two triangular 
graphs gives 

/«an,~4(^-a.j ^^^^ (24) 

and for the crossing box graph, similarly 

64 /4 y 

4rsbox ~ y asj (g2+^2+fc2)2- (25) 

These expressions have to be compared with the leading one-gluon term we have used (see Eq. (A8) of App. A) 



loGE ~ 167r^ "^'^ ■ ^^^^ 



Putting all things together, the overall error due to the omission of such NLO contributions to the BS kernel is then 



A/ _ fry -^triang ^ / -fcrsbox ^ 



and this produces 



AO AI 

~d T 



(28) 



By using Eqs. (20-23), after some algebra it is easy to recognize that the NLO effects on af'^P turn out to be of the 
same order, that is 

ti. AI 

Anlo^s ~ — , (29) 

which is what is assumed in the foregoing. The NLO errors do not exceed 5% for heavy quark states while they are 
enhanced up to 20% when light and strange quarks are involved. 

Finally, since the strength of the neglected coupling with other channels (OC) is obviously measured by the width 
of the state, one roughly estimates an error of the order of AttIq ^ ra/2 in the evaluation of nria ■ On this ground, 
for each determination of af^^ {QV) ^he related theoretical error is given by 

Tfl 

Usually the error Amexp on the experimental mass rriexp is much smaller than rtt/2 . When, however, this is not the 
case one has to consider also the experimental error AexpC^s , obtained from (30) by replacing rria Fa with 2mexp Amexp • 

Before discussing our results some comments are in order. 

First note that in the evaluation of Qa in (20) one has neglected the hyperfine splitting which however was taken 
into account in (23), bringing possibly to difl:erent values of 0;°^^ for the singlet and the triplet states (when there are 
reliable data for both). 
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Furthermore, the scmsitivity of the effective Q's determined from the specific coupHng (9), has been checked by 
analyzing the deviation of each Qo for a 25% shift of A^^^'^^g around the value 193 MeV, and one finds that the average 
change in the momentum scale amounts to 3%. This makes the resulting af^P{Q'^) reliable, at least qualitatively, 
even in the deep IR region (Q < 0.2 GeV), where the discrepancy with respect to massless a^^Q^) is sizable. 
There is a subtle point concerning the choice of the "unperturbed" ctg involved in the static Hamiltonian (17). 
Actually, the value adopted is very near to the a**^ pertaining to the bb{lS) state, but definitively smaller than the 
typical a^J^ 's. The point is that the hyperfine splitting is much more sensible than the c.o.g. mass to the behaviour 
of the unperturbed wave function at small distance (large momentum), which is specifically controlled by the value 
of the unperturbed as ■ As a result, the effective fixed value af^ in Eq. (17) that reproduces the same splitting as by 
using the coupling q:e(Q^) turns out to be significantly smaller than a^J^ calculated from the c.o.g. mass. Essentially, 
it was chosen a phenomenological value for the unperturbed as in order to have a good reproduction of the hyperfine 
splitting so as to reasonably reconstruct the c.o.g. of the doublet when one component is missing. It was then used 
the position of the c.o.g. (which is rather stable w.r.t. the unperturbed a^) to extract af^^{Q'^) . 
We finally stress again that the 1-loop analytic coupling with the above mentioned value of the scaling constant used in 
our computation, Eq. (9), differs by no more than 1% from the 3-loop analytic coupling in the region 0.5 < Q < 1.2 GeV 
where all the input states (tt, <p , J/ip and T ) fall. 



V. BS-MODEL RESULTS: CONCERT OF LOW AND HIGH ENERGY DATA VIA APT 

All results are displayed in details in tables I- VII of App. B, and pictorially in Fig. 3 taken from Ref. [1]. The first 
three columns specify the state and its experimental mass as given by [37]. The fourth column gives our theoretical 

results for the meson masses, and the last three give the effective Q , the relative 3-loop APT coupling a^\Q^) and 
the experimental coupling with errors (theoretical and experimental). 

In Fig. 3 values of a^^^ at the same Q from triplet and singlet states have been combined through a weighted average 
according to their errors (both experimental and theoretical). The c.o.g. of the light-heavy states (that are interpreted 
according to the j-j scheme) are also reported in the above figure. As can be seen, the agreement between the 3-loop 

analytic coupling with A^^^^g = 417 MeV, and the points representing our experimental values for the QCD coupling 
is quite good within the errors down to 200 MeV. 

At energies below 200 MeV a tendency of af^P{Q^) to diminish with Q seems to exist. As already noted, such a 
deep IR behaviour could be theoretically understood within theoretical models, in particular within the "massive" 
modification of APT in Sec. 2. Specifically, as displayed in Fig. 3, the one-loop couphng a{Q'^,m?) (13) with an 
efii'ective mass m^g ~ (38 ± 10) MeV reasonably fits all experimental points down to the very low Q region. 
Let us notice, however, that the analysis of such an extreme IR behaviour is based on high orbital excitations (D 
and F states), lying well above the strong decay thresholds and with large widths. As a consequence, the theoretical 
reliability of the method is lower at these scales, as apparent from the large estimated errors. Moreover, also the 
discrepancy between a^\Q'^) (used in the calculation) and a^\Q'^) (used as a reference term) rises above 10% at 
these scales. In fact, only the two states 772(1670) (interpreted as ss{l^D2)) and /2(2150) {ss{l^F2)), corresponding 
to Q ~ 120 MeV, generate af^'^ijQ'^) (marginally) out of the error bands, and the state /2(2150) (observed only once) 
has never been confirmed. 

Restricting our considerations to a sample of better established data, which excludes high orbital excitations as D and 
F states, the comparison with the BS meson masses yields a ~ 1 if an additional conventional error of 20 MeV is 
assigned to the latters. This error should account for the sources of theoretical uncertainty not explicitly evaluated, 

and produces an average additional error of roughly 5% on af'^P((5^) • 

At this point it is worthwhile to comment on the dependence of the results on the renormalization scheme. First 
of all, our definition of the coupling is implicitly contained in ansatz (1). Specifically, here one assumes both that 
ilnW^ is dominated by the OGE term after the subtraction of the area term, and that the OGE term is represented 
as the corresponding tree-level expression, the fixed coupling a^ being replaced with the running one a^{Q'^). The 
latter assumption amounts to the embodying all the dressing effects into ag((5^) (see, e.g.. Sec. 3.2 of Ref. [3] and 
references cited therein). It is worth noting that the coupling defined in such way is free of unphysical singularities 
by construction. At the same time, the analytic running coupling a-^{Q'^), which is involved in our calculations, is 
remarkably stable with respect to both the higher loop corrections and the choice of renormalization scheme (see 
Sec. 2 and detailed discussion of this issue in Ref. [30]). Thus one might expect that the same situation should also 
occur for Q!gxp(<3^), with the possible exception for the deep infrared region (see, e.g., Sec. 4.5 of Ref. [3] and references 
cited therein), where other nonperturbative effects could be relevant. 
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Notice that in our selection of states as a rule we have excluded irregular and incomplete multiplets. Of this type, e.g., 
in the light quark sector, are the 35' states {m-^ ?, — i is anomalously large and about twice as m2 ag^ — ni2 ) ' 
l^P {mi3pg being larger than mispj, l^D, F, G, H (incomplete). If however included in the analysis, all these 
states would bring the results in agreement with the general tendency outlined. 

Finally, Fig. 4 (taken from Ref. [1]) displays a synthesis of results for as{Q^) defined from bound states in the BS 
framework with high energy data. Here, low energy results are reported in a logarithmic scale from 100 MeV to 
220 GeV together with a sample of high energy data as given by S. Bethke [2], against the 3- loop analytic coupling 

a^\Q'^) (8) and its massive modification (13). Also shown in the figure is the common perturbative 3-loop coupling 
with IR singular behaviour that is definitively ruled out by the data. As can be seen, the BS-APT theoretical scheme 
allows a rather satisfactory correlated understanding of very high and rather low energy phenomena. 



VI. CONCLUSIVE REMARKS 

To summarize, we have exploited calculations within the Bethe-Salpetcr formalism adjusted for QCD, in order to 
extract an "experimental" strong coupling a'g^^{Q^) below 1 GeV by comparison with the meson spectrum. 
This work extends the analysis given in [1], providing technical details, the complete set of numerical results and their 
thorough discussion. 

A key point is the comparison of af^^ with the analytic coupling a^iQ'^) which avoids the hurdle of the unphysical 

singularities in the IR region [4] . 

The method consists in solving the eigenvalue equation for the squared mass operator as given by Eq. (2), obtained 
by a three dimensional reduction of the original BS equation. The relativistic potential U then follows from a proper 
ansatz (1) on the Wilson loop to encode confinement, and is the sum (2) of a one-gluon-exchange term Uqge and a 
confining term C/conf ■ The coupling occurring in the perturbative part of the potential needs to be IR finite since its 
argument has to be identified with the momentum transfer in the qq interaction, and this typically takes values down 
to few hundreds MeV. The usual perturbative running coupling has then been replaced by the 1-loop analytic 

expression ^ Eq. (9) with an effective QCD scale A^^^'^^^ = 193 MeV. This value reasonably reproduces the 3-loop 
analytic coupling, normalized at the Z boson mass along with world average [37] (i.e., A^"^-'^5 = 236 MeV which leads 
to A^^^^g = 417 MeV by continuous threshold matching, see App. C). 

Thus wc have taken advantage of our BS results for the meson spectrum, both in the light and heavy quark sector, 
to infer within this framework the fixed coupling value for each state that exactly matches the theoretical and 
experimental mass. 

Our results are twofold. On the one hand, as expected (given the good agreement of theoretical and experimental 
meson data), the 3-loop analytic coupling reasonably fits a'^'^^{Q^) from 1 GeV down to 200 MeV within the estimated 
theoretical and experimental errors, with very few exceptions. This confirms and yields a quantitative estimate of the 
relevance of the APT to IR phenomena down to 200 MeV. 

On the other hand, below this scale, the experimental points exhibit a tendency to fall under the APT curve. This could 
give a hint cither on the vanishing of q„(Q^) as Q ^ (in concert with some results from lattice simulations [17]), or 
on the existence of a finite IR limit of as{Q^) lower than the universal APT freezing value. The former IR behaviour 
can be theoretically understood in the framework of a recent "massive" modification [15] of the APT algorithm, 
which takes into account effects of a finite threshold in the dispersion relation. Since in the extremely low Q region 
confinement forces play the dominant role, the reasonable agreement between the "massive" APT model and the 
results of the BS formalism would suggest a relation between the linear potential, arising from the area term in the 
ansatz (1), and the thresholds effects in the analytic properties of the QCD coupling, to be further investigated. 
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APPENDIX A: SECOND ORDER BETHE-SALPETER FORMALISM 

In the QCD framework a second order four point quark-antiquark function and full quark propagator can be defined 
as 

Hi{xi,X2]yi,y2) = -^Trcoior(Ai(a;i, yi; A)A2(y2, 2:2; A)) (Al) 

and 

H2{x -y) = -^Trcoior(A(a;, y; A)) , (A2) 

{f[A]) = J DAMp[A] e^^«[^I/[A] , (A3) 

Mi? [A] = Det n^^;^[l + g^i^ An{i^^ djv — m^"")~^] and A(a;,y; A) is the second order quark propagator in an external 
gauge field. 

The quantity A is defined by the second order differential equation 

(D^Z?" + m2„„ - \g a^^F^,)A{x, y; A) = -d\x - y) , (A4) 

(a^'^ = |[7'',7''] and Dij_ = igA^) and it is related to the corresponding first order propagator by S{x,y;A) = 
{i'y^Dv + rucurr) A(x, y; A) , mcurr being the so-called current mass of the quark. 

The advantage of considering second order quantities is that the spin terms arc more clearly separated and it is 
possible to write for A a generalized Feynman-Schwinger representation, i.e., to solve Eq. (A4) in terms of a quark 
path integral [11,12]. Using the latter in (Al) or (A2) a similar representation can be obtained for H4 and H2. 

The interesting aspect of this final representation is that the gauge field appears in it only through a Wilson line 
correlator W. In the limit X2 ^ xi, y2 ^ yi or y ^ x the Wilson lines close in a single Wilson loop T and if T stays 
on a plane, ilnM^ can be written according to (1) as 

i\nW =^a^jdz'' j dz^'D^^iz - z') + (A5) 
ajdz'j dz^'Siz' - z°')\z - z'l J^' dx{l - [A^ + (1 - X)^?}' • 

The area term here is written as the algebraic sum of successive equal time strips and dz± = dz— {dz ■ r)r/r^ denotes 

the transversal component of dz. The basic assumption now is that in the center of mass frame (A5) remains a good 
approximation even in the general case, i.e., for non flat curves and when X2 ^ xi, y2 ^ yi or y ^ x. Then, by 
appropriate manipulations on the resulting expressions, an inhomogeneous Bethe-Salpeter equation for the 4-point 
function H/^{xi.X2',yi,y2) and a Dyson-Schwingcr equation for H2{x — y) can be derived in a kind of generalized 
ladder and rainbow approximation. This should appear plausible, even from the point of view of graph resummation, 
for the analogy between the perturbative and the confinement terms in (A5). 

In momentum representation, the corresponding homogeneous BS-equation becomes 

^ f dSi . / 1 k + u 1 ^ k + u 



X H^'> \^-P + kja'^^p{u)a'H'2 > [--P + kj , (A6) 

where a° = 1; a, 6 = 0, iii^; the cm. frame has to be understood, P = (tob, 0); $p(fc) denotes the appropriate second 
order wave function, that in terms of the second order field (j){x) = (17^!?^ + TOcurr)~^V'(2^) can be defined as the 
Fourier transform of (0|^(|)'^(— |)|P) . 
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Similarly, in terms of the irreducible self-energy, defined by H2{k) = i{k^ — Wcurr) ^ + ^(^^ ~ ™curr) ^ i^{k) H2{k) , 
the Dyson-Schwinger equation can be written 

m= J ^Jat{k-l;^,'^)a'^m)a» . (A7) 



The kernels are the same in the two Eqs. (A6) and (A7), consistently with the requirement of chiral symmetry limit 
[38], being given by 

Ufi{Q\P,p') = 167r^asPV^^a/3(Q) + 

+4a J d\e-^'^<\C\e{po)e{p'o) J\x{pIp'o' - [Xp'^PT + (1 - A^oP^^^ 



\C\Vp'o-pf° 



h;pa{Q;p,p') = -47ri^a«p«(^^Q<, - S^Qp)D^f3{Q) + 

w.MQ-.p.p') = ir5a.(«;o» - «e„)(*;<3<. - s;q,)6„>{Q} , (as) 

where in the second and in the third equation = has to be understood. Notice that, due to the privileged role 
given to the cm. frame, the terms proportional to a in (A8) formally are not covariant. 

In fact, it can be checked that r(A;) can be consistently assumed to be spin independent and Eq. (A7) can be rewritten 
in the simpler form 

V (277)4 p-m^ + niy ^^^^ 



with 

4 

{k-l)4k-lYD^^{k-l)-{k-iy'(k-IYD,,,Ak-l)] + 



R{k, I) = An-a, [{k + lY{k + l^D^^k - l)+ 



k_L and Ix denoting as above the transversal part of k and 1. Eq. (A9) can be solved by iteration resulting in an 
expression of the form r(A;^, k^), since (AlO) is not formally covariant. Then the constituent (pole) mass m is defined 
by the equation 

- m^u„ + r(m^ k^) = (All) 

and the dependence on k^ , being an artifact of the ansatz (1), is eliminated by extremizing m(k^) in k^ . 

The 3-dimensionaI reduction of Eq. (A6) can be obtained by a usual procedure of replacing H2{k) with z(fc^ — m^)~^ 
and lab with its so-called instantaneous approximation /^^^*(k, u) . In this way, one can explicitly integrate over uq 
and arrive to a 3-dimensional equation in the form of the eigenvalue equation for a squared mass operator Eq. (2), 
with [11] 



/.,lrrl.,/\ _ 1 / ^l +^2 finst,,, w'l + w'2 _a J> / 4 . 
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Finally by using Eq. (A8) one obtains Eqs. (14-16). 

Alternatively, in more usual terms, one could look for the eigenvalue of the mass operator or center of mass Hamiltonian 
HcM = M = Mo + V with V defined by M qV + VMp + = U. Neglecting term the linear form potential V can 
be obtained from U by the replacement / (^^+'^2){wi+w2) _^ — ^ — ^ ^j^^ resulting expression is particularly 

useful for a comparison with models based on potential. In particular, in the static limit V reduces to the Cornell 
potential 

V.t.t = ~-+cTr. (A13) 

Note that it is necessary to introduce a cut-off B in Eq. (A9). As a consequence the constituent mass turns out to 
be a function of the current mass and of -B, m = m (rricurr, B). Then if one uses a running current mass mcurr(Q^) 
we obtain a running constituent mass m{Q^) as it has been done in [13] (see also Ref. [24]). However the singular 
expression used there 

/ Q^/A^-l XTo/^/^o 

is not consistent with Eq. (9), and if a more consistent assumption is taken, e.g., 

„ / 1-11 o \7o/2/3o 

mcurr(Q') =m(4')(Q2-)j ^ (AI5) 

or the other resulting from the analytization of a similar expression with a^^\Q'^) replaced with the ordinary pertur- 

bative ai^\Q'^), the dependence of m on is strongly reduced. For this reason even the light quark mass is here 
treated as a constants to be adjusted with the the data** . 



**Iii this way the only role that is left to the DS equation is to justify the difference between the constituent and the current 
masses. 
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APPENDIX B: NUMERICAL RESULTS 

The tables below display the complete set of results as explained in Sec. 5. We recall the values of all the param- 
eters: a = 0.18 GeV^ K^^'^"^^ = 193 MeV, to, = 196 MeV [q = u,d) , = 352 MeV , = 1.516 GeV and 

mfc = 4.854 GeV. Meson masses are given in MeV. The last column displays the experimental coupling a'^^{Q'^^) 
with the theoretical error Anlo due to the next-to-leading order terms neglected, the theoretical error Ar from the 
half width r/2 and the experimental error Aexp respectively. 

t Center of gravity masses of the incomplete multiplets estimated in analogy with other multiplets. 

TABLE I. qq {q = u, d) 



States 


(MeV) 




mth 


Q 




± Anlo ± Ar ± Aexp 


I'So 




134.9766 ± 0.0006 1 
139.57018 ± 0.00035 J 


136 


401 


0.522 


0.534±0.122±-±- 


l^Si 


p(770) 


775.5 ± 0.4 


749 






0.517 ±0.122 d= 0.048 ±- 


lASS 




638 


613 








2^80 


7r(1300) 


1300 ± 100 


1223 


448 


0.502 


0.451 ± 0.114 ± 0.152 ± 0.081 


2^8i 


p(1450) 


1459 ± 11 


1363 






0.427 ± 0.114 ± 0.062 ± 0.010 


2A88 




159 


139 








l^Pi 


6i(1235) 


1229.5 ± 3.2 


1234 


209 


0.637 


0.688 ± 0.155 ± 0.124 ± 0.006 




7r2(1670) 


1672.4 ± 3.2 


1595 


144 


0.701 


0.544 ± 0.151 ± 0.364 ± 0.009 



TABLE 11. ss 



States 


(MeV) 


T/lexp 


mth 


Q 




af P ± Anlo ± Ar ± Aexp 


1^8i 


0(1020) 


1019.460 ± 0.019 


1019 


418 


0.514 


0.525 ± 0.098 ± 0.002 ± - 


2^81 


<^(1680) 


1680 ± 20 


1602 


454 


0.500 


0.435 ± 0.096 ± 0.068 ± 0.019 


l^Pi 


/n(1380) 


1386 ± 19 


1472 


216 


0.631 


0.824 ± 0.098 ± 0.083 ± 0.032 


1^P2 


/2(1525) 


1525 ± 5 "I 










l^Pl 


/i(1510) 


1518 ±5 \ 1521 


1484 






0.603 ± 0.098 ± 0.070 ± 0.009 


l^Po 


/o(1500) 


1507 ±5 J 










I'Dz 


r?2(1870) 


1842 ± 8 


1807 


149 


0.695 


0.658 ± 0.079 ± 0.318 ± 0.023 


1^F4 














l^Fs 




i 2165+ 


2070 


118 


0.733 


0.452 ± 0.064 ± 0.137 ± 0.024 


l^Fz 


/2(2150) 


2156 ±11 J 
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TABLE III. qs {q = u, d) 



States 


(MeV) 













exp ^ ^ ± Ar ± Aexn 


I'So 




497.648 ± 0.022 1 
493.677 ± 0.016 J 


491 


409 


0.518 


0.529 ±0.122±-±- 


l^Si 


( if* (892)° 
1 K*(892)^ 


896.00 ± 0.25 
891.66 ± 0.26 


1 893.11 


887 






0.526 ±0.122 ±0.017 ±- 


lASS 




398 




396 








2^81 


A" (1410) 


1414 ± 15 


1485 


451 


0.501 


0.571 ± 0.117 ± 0.102 ± 0.013 


l^Pi 


ii:i(1270) 


1272 ± 7 


1355 


213 


0.634 


0.820 ± 0.129 ± 0.081 ± 0.012 


1^P9 f 
1 


1432.4 ± 1.3 1 












l^Pl 


A'.p (1430)=*= 
" ifi(1400) 


1425.6 ± 1.5 
1402 ± 7 1 


1 1417.7 


1367 






0.583 ± 0.129 ± 0.133 ± 0.007 


1 Po 


7^0* (1430) 


1414 ± 6 J 












l^Da 


isTKlTSO) 


1776 ±7 ] 














/!'2(1770) 


1773 ±8 


> 1763 


1712 


150 


0.694 


0.617 ± 0.113 ± 0.273 ± 0.031 


l^Di 


/i'*(1680) 


1717 ±27 J 












1^F4 


/s:4*(2045) 


2045 ±9 ] 














7^3(2320) 


2324 ±24 


> 2121 


1973 


116 


0.736 


0.248 ± 0.095 ± 0.413 ± 0.071 


l^Fa 


is:2*(1980) 


1973 ±25 J 
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TABLE IV. cc 



States 


(MeV) 






mth 


Q 




af P ± Anlo ± Ar ± Aexp 


I'So 
l^Si 
lASS 
2^80 
23Si 
2ASS 
3^81 
4^Si 


?7c(15) 

V(2S) 

V'(4040) 
i/)(4415) 


2980.4 ± 1.2 
3096.916 ± 0.011 
117 
3638 ± 4 
3686.093 ± 0.034 
48 

4039 ± 1 
4421 ± 4 


2980 
3097 

118 
3595 
3653 

58 
4030 
4337 


561 

500 

483 
474 


0.464 

0.483 

0.489 
0.492 


0.467 ± 0.025 ± 0.008 ± 0.001 
0.467 ± 0.025 ± - ± - 

0.446 ± 0.023 ± 0.007 ± 0.004 
0.455 ± 0.023 ± - ± - 

0.485 ± 0.022 ± 0.049 ± 0.001 
0.384 ± 0.022 ± 0.042 ± 0.006 


I'Pi 

l^Pl 

l^Po 
2 3P2 

2 3Po 


hc{lP) 

Xc2(lP) 
Xcl(lP) 

Xco(lP) 

Xc2(2P) 

X(3872) 


3525.93 ± 0.27 
3556.20 ±0.09 \ 
3510.66 ±0.07 \ 3525.3 
3414.76 ±0.35 J 
3929 ±5 1 

\ 3915+ 
3871.2 ±0.5 J 


3532 
3537 

3929 


269 
274 


0.592 
0.589 


0.631 ± 0.012 ± - ± - 
0.640 ± 0.012 ± 0.002 ± - 

0.644 ± 0.013 ± 0.027 ± 0.009 


l^Di 
2^D3 
2^D2 
2^01 


V'(3770) 


3771.1 ±2.4 J 


> 3820+ 


3822 


190 


0.654 


0.707 ± 0.008 ± 0.030 ± 0.006 


V'(4160) 


i 4183+ 
4153 ±3 J 


4150 


198 


0.646 


0.606 ± 0.009 ± 0.132 ± 0.008 



Quantum numbers of hc{lP) and X(3872) mesons are not well established. 



TABLE V. bb 



States 


(MeV) 


'^exp 




mth 


Q 




aT^ ± Anlo ± Ar ± Aexp 


l-'Si 


T(15) 


9460.30 ± 0.26 




9461 


951 


0.381 


0.378 ± 0.006 ± - ± - 


2^81 


T(2S) 


10023.26 ± 0.31 


9987 


630 


0.445 


0.416 ± 0.004 ± - ± - 


3*Si 


T(3S') 


10355.2 ± 0.5 




10321 


552 


0.466 


0.433 ± 0.003 ± - ± - 


43Si 


T(45) 


10579.4 ± 1.2 




10588 


517 


0.478 


0.493 ± 0.003 ± 0.013 ± 0.002 




T(10860) 


10865 ± 8 




10820 


497 


0.484 


0.424 ± 0.003 ± 0.078 ± 0.011 




T(11020) 


11019 ± 8 




11034 


506 


0.481 


0.508 ± 0.003 ± 0.057 ± 0.012 


1^P2 


X62(1P) 


9912.21 ±0.26 ±0.31 ] 












l^Pi 


Xbi{lP) 


9892.78 ± 0.26 ± 0.31 


1 9899.87 


9880 


387 


0.528 


0.519±0.002±-±- 


l^Po 


Xm(1P) 


9859.44 ± 0.42 ± 0.31 J 














X62(2P) 


10268.65 ± 0.22 ± 0.50 1 














X6i(2P) 


10255.46 ± 0.22 ± 0.50 


1 10260.24 


10231 


343 


O..549 


0.524 ± 0.002 ± - ± 0.001 


2"Po 




i()2:!2.r, ± 0.1 ± ()..-,() J 
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TABLE VI. Light-heavy quarkonium systems 



States 


(MeV) 








mth 


Q 




Qf P ± Anlo ± Ar ± Aexp 


qc 


















I'So 


\ D" 


1869.3 ± 0.4 1 
1864.5 ± 0.4 J 


|- 1867.7 


1843 


459 


0.498 


0.488 ± 0.082 ± - ± - 


1 3q J 


\ D*(2010)± 
[ D*(2007)° 


2010.0 ± 0.4 1 
2006.7 ± 0.4 J 


\ 2008.9 


zUUU 






u.4yy ± u.uttz ± — ± — 


lASS 




141 ± 1 




157 








IP2 j 

- i 


f D2*(2460)± 
[ D|(2460)° 
f Di(2420)± 
[ i?i(2420)° 


2459 ± 4 ] 
2461.1 ±1.6 1 
2423.4 ± 3.1 1 
2422.3 ±1.3 J 


1 


2443 


232 


0.619 


0.651 ± 0.051 ± 0.028 ± 0.005 



qh 

l^So 5279 4 ± 5 } ^^^^'-^ ^^^^ ^'^"^^ 0.456 ± 0.036 ± - ± - 

l^Si B* 5325.0 ± 0.6 5311 0.471 ± 0.036 ± - ± - 
lASS 46 ± 3 64 



TABLE VII. Light-heavy quarkonium systems 



States 


(MeV) 




JTlth 


Q 




QfP ± Anlo ± Ar ± A^xp 


sc 
l^So 




1968.2 ± 0.5 


1959 


472 


0.493 


0.494 ± 0.055 ± - ± - 


l^Si 


Dt 


2112.0 ± 0.6 


2109 






0.497 ± 0.055 ± - ± - 


lASS 




144 ± 1 


149 








IP2 
IPi 


73,2(2573)* 
D,i(2536)* 


2573.5 ±1.7 1 
2535.35 ± 0.34 ± 0.5 J 


2545 


236 


0.616 


0.626 ± 0.036 ± 0.009 ± 0.002 


sh 
l^So 




5367.5 ± 1.8 


5343 


535 


0.472 


0.457±0.024 ±- ±0.001 


l^Si 




5412.8 ± 1.7 


5408 






0.473 ± 0.024 ± - ± 0.001 


lASS 




17 ± 1 


05 








IP 


B:^(5850) 


5853 ± 15 


5830 


255 


0.602 


0.592 ± 0.013 ± 0.042 ± 0.027 
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APPENDIX C: 3-LOOP ANALYTIC COUPLING 



As mentioned, we compare our results with the 3-loop analytic coupling (8), the difference with respect to the 4-loop 
approximation being negligible [30]. The exact 3-loop spectral density (7) has a rather cumbersome structure in terms 

of the Lambert function (sec Rcf. [33]). However, for practical purposes, it can be approximated with a high accuracy 
by the discontinuity of the 3-loop perturbative running coupling (as given by PDG [37]) across the physical cut [25]: 

= ^ ^2+^2)3 W^' - - (3*' - ^')m] ■ (CI) 



In this equation t = ln{a/A^), 

J{t) =2t- Bi[tGi{t) + G2{t)] + BfGi{t)[2G2{t) - 1], (C2) 

R{t) =e-TT^ - Bl[tG2{t) - TT^Glit)] + 

BfiGlit) - TT^Glit) - G2{t) - 1] + B2, (C3) 



Gi{t) = ^ - - arctanf-) , G2(t) = ^ Ht^ + n^) (C4) 

with Bj = /3j//3o^^ being the combination of the /3 function expansion coefficients 

Po = ^ f 11 - In!) , (C5) 



In Fig. 5 the 3-loop analytic coupling (8) is compared to the same level perturbative expression, both normalized 
at the Z boson mass to the world average value as(Mf ) = 0.1176 ± 0.0020 [37] and evolved at the heavy quark 
thresholds crossing by continuous matching conditions. This gives for the analytic coupling the scaling constant 
A^^Lg ~ (417 ± 42) MeV in the IR domain. 
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FIG. 1. Quarkonium spectrum, three different calculations. Diamonds refer to the freezing prescription for the running 
coupling, squares and circles refer to the calculation with the 1-Ioop analytic coupling (9) and two different expressions for 
running constituent masses of light quarks, a solution of the Dyson-Schwinger equation and a phenomenological function of the 
cm. quark momentum respectively. Horizontal lines represent experimental data. 
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FIG. 3. Comparison between points extracted from Tables I- VII in App. B and the three-loop analytic coupling (8) with 
A^^' 3 = (417 ± 42) MeV (solid curves). The "massive" one- loop analytic running coupling (13) (dashed curves) refers to 
^(i,eff) ^ 204 MeV and the effective mass m^ff = (38 ± 10) MeV. Three-loop perturbative coupling (dot-dashed curve) corre- 
sponds to A^^^'^3 = 318 MeV. Circles, stars and squares refer respectively to qq , ss and qs with q = u,d (light-light states). 
Heavy-heavy states, cc and bb , are represented by diamonds and crosses. Finally in the light-heavy sector, asterisks stay for 
qc and qb , whereas plus signs stand for sc and sb . Data for triplet and singlet states referring to the same multiplet are 
combined in a weighted average according to their errors. Error bars include both theoretical and experimental uncertainty 
and are drawn only if relevant. 




FIG. 4. Summary of low (o) and high energy (o) data against the three-loop analytic coupling (8) (solid curve) and its 
perturbative counterpart (dot-dashed curve) both normalized at the Z boson mass. Also shown is the "massive" one-loop 
analytic coupling (13) (dashed curve) same as in Fig. 3. 
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FIG. 5. Three-loop QCD coupling below 2 GeV: the analytic expression (8) (solid curve) versus the perturbative one 
(dot-dashed curve), both normalized at the Z boson mass according to the world average [31]. 



